Evidence for a charge Kondo effect in Pbi.^^TlajTe from measurements of 

thermoelectric power 
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We report measurements of the thermoelectric power (TEP) for a series of Phi-xT^xTe crystals 
with X = 0.0 to 1.3%. Although the TEP is very large for x = 0.0, using a single band analysis 
based on older work for dilute magnetic alloys we do find evidence for a Kondo contribution of 11 - 
18 i-iV/K. This analysis suggests that Tk is ~ 50 - 70 K, a factor 10 higher than previously thought. 

PACS numbers: 72.15.Qm, 72.15.Jf, 74.70.Ad 



I 

o 
o 



o 



The traditional "spin" Kondo effect occurs in a non- 
magnetic host metal containing a small concentration 
of magnetic impurities. Here the anti-ferromagnetic 
exchange interaction (J) between the local magnetic 
moment and the conduction electrons gives rise to a 
— I logT term in the temperature (T) dependent elec- 
trical resistivity p{T) and other unusual propertiesi^i^ 
including an anomalously large and T-dependcnt ther- 
moelectric power (TEP or S). Below the Kondo temper- 
ature (Tk) which can be extremely small, the spin of the 
impurity is compensated by a cloud of conduction elec- 
tron spins extending over a distance ~ hvp/kBTK, where 
Ui? is the Fermi velocity of the host metal. There is a 
"triple peak" in the impurity density of states (DOS)"^, 
with two side lobes derived from the spin-split virtual 
bound state and a narrow peak of width ksTK at the 
Fermi energy. It is thought^ that a similar description 
applies to the charge Kondo effect discussed here. 

Experimental^ and theoretical^ evidence has been re- 
ported for a charge Kondo effect in Pbi_2;Tla;Te crys- 
tals with X > 0.3%. It arises because T1+ and T\^+ 
ions both have filled shells, i.e. 6s^ and 6s°, that can 
be more stable than the Gs^ state of TP+. So when Tl 
is in an environment favorable for divalency there can 
be two degenerate, or nearly degenerate, non-magnetic 
charge states Tl"^ and TP+ that can be described by the 
Anderson model and can also give rise to a Kondo effect^. 
The resulting charge fluctuations are thought to be im- 
portant for superconductivitji^i^i^ in Pbi_2;Tlj;Te. In the 
present paper we report TEP data for single crystals of 
Pbi_a;Tla;Te and discuss evidence for an anomalous con- 
tribution associated with a charge Kondo effect. 

The crystals of Pbi_^TUTe with x =0, 0.2, 0.3, 0.6, 
1.1 and 1.3% were from the same preparation batches 
studied previously^'^'^. Their thermopower was mea- 
sured using CuBe wires^, whose TEP was measured sep- 
arately relative to a cuprate superconductor and found 



to be small with a maximum of 0.2 iiV/ K near 30 K. Ini- 
tially the temperature gradient, AT, was measured using 
a constantan-chromel thermocouple made of 25 fj,m di- 
ameter wires, and glued to the mm-sized crystals with 
GE varnish, but for the data reported here very small 
diode thermometers^ were used. These were less straight- 
forward to mount, but gave more reliable measurements 
of AT, especially at low T where the sensitivity of the 
thermocouple decreases. The diodes were also attached 
to the sample using GE varnish, and the distance be- 
tween them, measured to ±3 to 5 % using a binocular 
microscope, is the main source of error in the TEP. 

For many dilute metallic alloys p(T) approximately 
obeys Matthiessen's rule which states that the contri- 
butions pj arising from two or more different scattering 
mechanisms (j = 0,1,2..) are simply additive. It gener- 
ally holds reasonably well when T is comparable to the 
Debye temperature (Od) or when the impurity resistivity 
is smaller than that caused by electron-phonon scatter- 
ing, but at lower T there can be significant deviations. 
One probable reason is that conservation of crystal mo- 
mentum (k) is less strict in an alloy because the finite 
electron mean free path limits the size of electron wave 
packets^. The equivalent Nordheim-Gorter (NG) rulei^ 
for combining contributions Sj to the total electron dif- 
fusion TEP S"^ from different scattering mechanisms in 
the same band is: 
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It relies on Matthiessen's rule being obeyed and on 
the scattering being effectively elastic, which is only true 
for T ^ Qo or in the residual resistivity (pres) region. 
Resistivity data^- for Pbi_j;Tlj;Te crystals from the same 
preparation batches are shown in Fig. 1. Matthiessen's 
rule is obeyed reasonably well for x > 0.6% for all T and 
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for X = 0.3% for T < 150 K. The insert shows that p^es 
is small for x < 0.2% but then increases linearly with x 
for X > 0.3%. This is consistent with other evidence that 
valency skipping sets in near x = 0.3%^ and, together 
with Matthiessen's rule, implies that there are no gross 
changes in electronic structure at higher x. The linearity 
also suggests that the Tl impurities act as independent 
scattering centers for x > 0.3%. TEP data measured 
for the six samples are shown in Fig. 2. All Tl-doped 
samples have S ~ 100 - 140 ^V/ K at 300 K with similar 
slopes and a positive curvature at lower T, while the pure 
PbTe sample has a larger TEP, 300 pV/K at 300 K, and 
a negative curvature. 




T(K) 

FIG. 1: Color online: resistivity vs. temperature for 
Pbi-2,Tlj;Te crystals^ from the same preparation batches as 
those studied here, values of x given in %. The insert shows 
residual resistivity vs. x for all crystals measured, p(T) data 
for all X values in the insert are given in Ref 13. 



Analysis of various Kondo alloys (e.g. AuFe, AuMn 
and AuCr)iiii^ suggests that S"' in Eqn. 1 contains three 
terms, 5*0, 5*1 and Sk- Sq arises from the energy depen- 
dence of the potential scattering term V in the Kondo 
Hamiltonian, and is related to the energies and widths 
of the 3d virtual bound states. It is linear in T and in 
Eqn. 1 is weighted by Pres = Po- The electron-phonon 
scattering term, 5*1, is also expected to be linear in T 
both for T > Gd/S and T < Gn/lO but with a factor 
of 3 smaller slope at low T. For PbTe, heat capacity 
data^ give Qd =168 K. So , Si, po and pi are often cal- 
culated using first order perturbation theory, while Sk, 
the Kondo TEP contribution, arises from higher order 
scattering processes involving non-cancelling Fermi fac- 
tors. It is also weighted by po in Eqn. 1. A broad peak 
in 5''' is often observed near the Kondo temperature Tr-, 
and often used to estimate Ti^^. However it has been 



argued that Sk could be constant for T > TKr^^, as 
predicted by high T perturbation theory^ and that the 
fall above the peak is caused by the other T-dependent 
terms in Eqn. 1. For T < 0.1 - 0.15Tk^, Sk falls to 
zero as T\ and fits Sk = AT/{T + 0.35Tk) for T < Tk. 
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FIG. 2: Color online: measured values of the TEP for 
Pbi_2,TliTe crystals with x = 0.0 and 0.2% -left hand scale 
and X — 0.3, 0.6, 1.1 and 1.3 % - right hand scale. The inset 
shows the weighted Kondo contribution, Skw = SkPo/p{T) 
for the values of x{%) shown. 

In view of the above discussion, we expect Sk to be 
constant or to fall slowly with T above a cut-off temper- 
ature To > Tk, and so above Tq Eqn. 1 gives: 

^ ^ aTpo + (3Tpi + 7po 
po + pi{T) 

where a, (3 and 7 are independent of x and T. Fits to 
Eqn. 2 were made from 300 K to Tq = 90, 70, 50 or 30 
K, with a, (3 and 7 as free parameters and po and pi{T) 
obtained by fitting the data in Fig. 1 to p ~ A + BT + 
CT^. Values of a, /3 and 7 and A, B and C arc given 
in Table 1. Tq =70 K gave shghtly better fits. Fits from 
150 to 70 K gave similar results and for x = 0.3% we 
used this range because of the curvature in p{T) shown 
in Fig.l for this sample. Below 70 K, Sk{T) was obtained 
from Eqn. 2 by using the formula Sk{T) = Smeas p/ Po ~ 
aT — (3Tpi I Po and is shown in the main part of Fig. 3a. 
Our fitting procedure should give a constant value for 
Sk above 70 K. The variations seen in Fig. 3a arise from 
residual errors in the fits multiplied up by p/po and are 
only a few times the noise level. 

For an isotropic parabolic band with Fermi energy Ep 
P ~ TT^k%/{eEF) for T > Qd/^ and n^kl/SEp for 
T < ez5/10i°'i&, i.e. T < 17 K for PbTe. The /3 values in 



3 



TABLE I; Fitting parameters, A, B, C, for p{T) and a, {3, 7, 
for TEP. 





A 
mClcm 


B 

mflcm/K 


C 

mflcm/K^ 






7 

liV/K 


0.3 


0.2017 


1.080x10"^ 


9.29x10-" 


-0.581 


0.446 


14.3 


0.6 


0.4124 


2.315x10"^ 


3.37x10-" 


-0.362 


0.682 


11.5 


1.1 


0.7371 


3.343x10"^ 


5.97x10-' 


-0.211 


0.837 


12.1 


1.3 


0.8472 


3.640x10"^ 


-2.6x10-** 


-0.211 


0.947 


17.7 



Table I give Ep between 73 to 160 raeV so their sign and 
magnitude are reasonably consistent with Ep being mea- 
sured relative to the bottom of the S band^, despite the 
fact that the hole pockets are actually ellipsoidal. But 
the changes of (3 with x are not understood since Hall 
data^ suggest a slight increase in hole concentration, i. 
e. in Ep with x. The magnitudes of a are similar to 
(3 but of opposite sign corresponding to higher energy 
holes being more strongly scattered by the Tl impurities. 
Within a Kondo picture this arises from the asymmetry 
of the Friedel-Anderson virtual bound states describing 
the valency fluctuations, but a low value of Ep for the S 
band could also affect a. Finally the values of 7 ranging 
from 11 to 18 jiV/K are similar to those associated with 
the traditional Kondo effect and the asymmetric Ander- 
son model. 

Fits to Eqn. 2 for the a; = 0.0 and 0.2 % samples 
were not good and are not shown. As shown in Fig. 1 
for these two samples pres values are very low and p{T) 
much more curved than for x > 0.3%. These differences 
and Hall dats*^ suggest that the L band dominates elec- 
tronic transport for low x while for x > 0.3% the E band 
plays a leading lole^. The TEP of the 0.2 % crystal 
varies as AT between 40 and 100 K with A = 0.2 nV/K"^ 
but rises sharply at higher T. This value of A is rea- 
sonably compatible with Ep = 200 meV obtained from 
analysis of Hall data^ but the strong increase at higher T 
is not understood. The deviations from linearity for x — 

0. 2 % shown in Fig. 3b are ascribed to the gradual onset 
of phonon drag below 40 K. This is considerably lower 
than for most metals where it is usually ~ Qpr^ and 
may arise because both the L and S hole pockets have 
small Fermi wave- vectors (kp). So electron-phonon scat- 
tering within a pocket will be suppressed to lower T until 
typical phonon wave-vectors ~ (T/0£))(7r/a), where a is 
the lattice spacing, become smaller than 2kp. We can- 
not assume that phonon drag corrections for the S band, 

1. e. X > 0.3%, will be similar to those for x = 0.2%. 
Perhaps the only way to obtain the phonon drag term 
for X > 0.3% is to suppress superconductivity by apply- 
ing a magnetic field and measure to much lower T where 
it will eventually vary as T^. Below T ~ 8 K the raw 
data in Fig. 3b show an approximately T-linear, and 
x-independent decrease that is typical of the Kondo ef- 
fect for T < 0.1 - 0.2 Tp^. However according to the 
preceding analysis the aT terms shown by the dashed 
lines in Fig. 3b are still present at low T. They should 
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FIG. 3: Color online: (a) Kondo contribution to the TEP 
of Pbi-ajTliiTe single crystals with values of x shown in %. 
The insert shows T /Sk vs. T plots used to find Tk from the 
empirical, low T, lawii, Sk = AT / [T + Q.Z^Tk)- (b) Raw 
TEP data for the values of a; in % shown. A solid line shows 
S'^ for X = 0.2% and dashed fines show aT for x > 0.3% 



be subtracted to obtain Sk(T) and this spoils the x- 
independence to some extent. 

We estimated Tk using two methods. The first one is 
based on the widely held view^iii^ that that there is an a;- 
independent peak in the TEP at Tk in dilute magnetic 
alloys, where the host metal has a very small TEP. If 
the TEP of the PbTe host metal and Sq were both small 
then we would measure a weighted value Skw = SkPo / P- 
Plots of Skw found in this way, using the Sk{T) data 
in Fig. 3 are shown in the insert to Fig. 2. We see that 
there are indeed broad peaks near 60 K for x = 0.6, 1.1 
and 1.3%. The second method is based on the empirical 
law 5"^ = AT/{T + 0.35Tk>M for which plots of T/Sk 
vs. T give a straight line extrapolating to Sk = at 
-0.35Tk- As shown in the insert to Fig. 3a such plots are 
reasonably linear and give Tk values ranging from 45 to 
75 K for the 4 samples. The deviations below 20 K for 
X = 1.1 and 1.3% could arise from phonon drag. Values 
of Tk obtained in these two ways are significantly higher 
than Tk ^ 6 K^ which was estimated by fitting resistivity 
data to p = p(0)[l - {T/Tk)^] but with considerable 
uncertainty in the appropriate value of p{0) . On the other 
hand, using the fact that the T^ law normally extends up 
to 0.1 Tkt^: gives Tk ^ 40 K - in closer agreement with 
our estimates from the TEP. 
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Although we should be cautious about applying spin 
Kondo formulae to a doped semiconductor, PbTe, the ob- 
servation of Matthiessen's rule is consistent with a single- 
band metallic picture and minor changes in hole concen- 
tration above x — 0.3%. Taking Tk — 60 K^^ does affect 
some of the previous conclusions^. For example on the 
basis of the measured specific heat coefficient 7(0) and 
the depth of the resistivity minimum it was suggested^'^ 
that only x^ff/x ^ 0.01 of the Tl impurities were degen- 
erate to within Tk = 6 K and hence only these were effec- 
tive Kondo scatterers. With Tk = 60 K, Xeff/x becomes 
~ 0.1 and also the radius of the charge cloud is smaller, 
~ 25 nm, taking kp = 10^ cm~^ and m = O.Gwe, i-e. vp 
= 2 lO*" cm/sec. Accordingly at x ^ 1%, on average there 
would still be about 10 other Tl atoms within the charge 
cloud of one impurity, but because Xeff/x 0.1 only 
one of these would be a Kondo scatterer. We note that 
the charge clouds must overlap and be delocalized at low 
T (as are the Kondo states in heavy Fermion supercon- 
ductors) because otherwise they would not contribute to 
superconductivity and there would be a residual 7(0) in 
the superconducting state. Such considerations of charge 
overlap may make detailed analysis in terms of the L and 
S bands of undoped PbTe less straightforward than pre- 
viously thought. 

Assuming that there is also a narrow peak ^ ksTK 
wide in the DOS for charge Kondo impurities, then 
the field scale for magneto-transport effects is larger 
^ ksTK/ fJ-B- Also one would expect a magnetic suscep- 
tibility contribution xi'T) = IJ^^Xef f N av I ksiT + Tk) 
emu/mole where Nav is Avogadro's number, /i^, the 
Bohr magneton, and from above a^e// ~ 0.1a;. This would 



be undetectable in the available xC^) data^ii but might 
show up in detailed measurements and analysis, like those 
made for AlMn alloya^^. There are some indications that 
Xeff/x could be larger than 0.1. Using the formula^ 
7(0) = 0.41287r2fcB/(6Tif) gives 7(0) = 0.94 mj/mole/K^ 
for X — 1% compared with the average slope from the 
raw experimental dataA of 7(0)/x = 0.57 mj/mole/K^/% 
implying that Xeff/x ~ 0.6. We also note that the 
measured values of Ax(0)/A7(0) are remarkably close 
to 6/i^/7r^fc^, the value expected from spin Kondo the- 
ory, namely a Wilson ratio of 2 rather than 1 for non- 
interacting fermions. However for Xeff/x ~ 0.6 it is 
difficult to understand the T dependence since the x(T) 
curves for x between 0.3 and 1.3 % are essentially parallel. 
Although the effective mass and band gap of PbTe are 
T dependent-lS , it seems unlikely that this would exactly 
compensate any 1/(T + Tk) behavior from the Kondo 
effect. 

In summary, the large TEP of the host material makes 
the data analysis less straightforward than for noble 
metal-based Kondo alloys. Using a provisional single- 
band picture we do see reasonably clear evidence for a 
Kondo effect but with Tk a factor of 10 or so larger than 
the earlier work. Band calculations of the TEP of PbTe 
as a function of hole doping might provide a further test 
of this conclusion. 
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